function LMP=marginalMetLaplace(H,logPosterior)
H=0.5*(H'+H);
[~,flagNotPD]=chol(H);
if flagNotPD~=0
    LMP=-inf;
    return
end
[~,vv,~]=svd(H);
logDet=sum(log(diag(vv)));
NP=size(H,1); 
LMP=0.5*NP*log(2*pi)+0.5*logDet+logPosterior; 
sprintf('Marginal Likelihood=%10.5f',LMP); 